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High-T QCD and dimensional reduction: measuring the Debye mass* 
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We study the high-temperature phase of SU(2) and SU(3) QCD using lattice simulations of an effective 3- 
dimensional SU(A'^) + adjoint Higgs -theory, obtained through dimensional reduction. We investigate the phase 
diagram of the 3D theory, and find that the high-T QCD phase corresponds to the metastable symmetric phase 
of the 3D theory. We measure the Debye screening mass mo with gauge invariant operators; in particular we 
determine the 0{g^) and 0(17^) corrections to mn- The corrections are seen to be large, modifying the standard 
power-counting hierarchy in high temperature QCD. 



The Debye mass (inverse screening length of 
color electric fields) characterizes the coherent 
static interactions in QCD plasma, and its numer- 
ical value is essential for phenomenological dis- 
cussion of the physics of the QCD plasma. For 
SU(iV) QCD, iV = 2, 3, with Nf massless quarks, 
the Debye mass can be expanded at high temper- 
atures in a power series in the coupling g: 



turbative analysis. It retains the essential in- 
frared physics of the original theory, and since it 
is bosonic even for Nf > 0, it can be studied eco- 
nomically with lattice Monte Carlo simulations. 
Recently it has been very successfully applied to 
the Electroweak phase transition |7j . 

The effective action is derived with the Green's 
function matching technique 
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where the leading order perturbative result is 
TO^O = [N/?,+Nf/&Y/'^gT. The logarithmic part 
of the 0{g^) correction can be extracted pertur- 
batively [Q, but cat and the higher order correc- 
tions are non-perturbative. Our aim is to evalu- 
ate numerically the coefficients cat and dN,Nf A 
detailed discussion of the results can be found in 

Hi. 

An effective 3D action, obtained through di- 
mensional reduction is a powerful tool for 
studying high-T QCD. The effective action can be 
derived perturbatively without the infrared prob- 
lems associated with the standard high-T per- 



^0 is a remnant of the temporal gauge fields and 
belongs to the adjoint representation of SU(7V). 
The Lagrangian (g) gives the Green's functions 
to a relative accuracy 0{g'^) sufficient for the 
accuracy of the expansion in cq. (|^). 

For Nf = 0, the 3D couphngs gl, 7n\ and \a 
are related to the temperature T (and the 4D 
gauge coupling g^{p), which is evaluated at the 
optimized scale ^ « 27rT) by 
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The dynamics of the 3D theory is fully character- 
ized by the dimensionless ratios x and y above 
and the dimensionful gauge coupling g^. The 
presence of fermions only modifies the numerical 
factors in eqs. (|^-||) . 

Due to the superrenormalizability of the 3D ac- 
tion the continuum<->lattice relations of the cou- 
plings become very transparent (for a detailed 
discussion, see In particular, the lattice 

gauge coupling is related to the continuum 
gauge coupling and the lattice spacing a by 
l3G = 2N/{gla). 

The phase diagram of the SU(2) -I- adjoint 
Higgs theory is shown in Fig.|l| (for convenience, 
plotted in the {x , xy)-plane) . For SU(3) the di- 
agram is qualitatively similar. At small x, the 
transition is very strongly first order, but becomes 
rapidly weaker when x increases. At a; w 0.32 
there is a critical point, after which only a cross- 
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Figure 1. The phase diagram of the 3D SU(2) -|- 
adjoint Higgs theory. The plot symbols show the 
location of the transition given by the simulations 
(extrapolated to V ^ oo), and the y = yc{x) 
curve is a 4th order polynomial interpolation of 
the data. Dashed line indicates where the transi- 
tion becomes a cross-over. The straight lines are 
the 'dimensional reduction lines' ydr(a^) for the 
number of fermion flavors indicated. 



over remains. The Nf = ydr(a^) -line is given 
in eq. ^ . The 3D theory is well defined on the 
entire (x, ?/)-plane, but only along ydrix) does the 
3D theory describe the physics of the 4D SU(2) 
gauge theory. Along this line x is related to the 
temperature as shown on the top axis of the plot 
(for SU(2), A5^«1.2r,). 

In the physically relevant region T ^ Tc ~ 
Aj^ the ydr line is in the broken phase . How- 
ever, the 3D theory cannot describe 4D high-T 
physics in the broken phase, since the pertur- 
bative 4D<->3D connection is not valid there 
Thus the symmetric phase is the physical one. 
Due to the strong 1st order nature of the transi- 
tion at small x, the symmetric phase is strongly 
metastable (shown with the dashed line in Fig. |^). 
If initially prepared to be in the symmetric phase, 
the system remains there for the duration of any 
realistic Monte Carlo simulation. 

The Debye mass can be defined as the mass 
of the lightest 3d state odd under the reflection 
^0 —^0 The lowest-dimensional gauge 
invariant operator fulfilling this is 



hi = EijkTrAoFjk 



(7) 



To enhance the signal, we measure the correla- 
tion function {hi{0)hi{z)) using several recursive 
blocking levels . We perform the measurements 
along the metastable ydi{x) -lines, eqs. (^,^. The 
results are shown in Fig. ^, in units of 4D gT 
(= gl^-iy/N in 3D units). The Monte Carlo 
runs are performed with several lattice spacings 
a (Pg)- The top scales in Fig. ^ show the physi- 
cal temperature T/Aj^ along ydri^) -lines. Note 
that the highest temperatures are larger than 
IQioo xAj^^ IQi™ X T,. 

We fit the data to the 2-parameter ansatz in 
eq. (|^). At small x (large y), the quality of the 
fits is very good. We use only data in the range 
X < 0.08, so that the horizontal plateaus at large 
X are excluded from the fits. The results of the 
fits are l3l 



SU(2): 
SU(3): 



C2 = 1.58(20) 
C3 = 2.46(15) 



62 = -0.03(25) 

63 = -0.49(15), 



(8) 



where bw = dN,Nf\/N/3 + Nf/6. For iV = 2 we 
can only verify that d2.Nf is close to zero. Note 
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Figure 2. The Debye mass for SU(2) (top) and 
SU(3) (bottom), as functions of a; (top scales show 
T/Aj^ through eq. @)). The dashed hne indi- 
cates the leading order m^°, and the continuous 
line the 2-parameter fit to eq. (|l|) with the pa- 
rameters as in eq. m). 



that writing cjv = Ncn, one has cjv = 0.79 ±0.10 
(iV = 2), 0.82 ±0.05 (A^==3). 

The leading contribution to mo is dominant 
only at extremely large T - indeed, for SU(3) 
the leading term is larger than the 0{g^) correc- 



tion for r/A^> 10^, implying that the leading 
term only dominates when QCD anyway merges 
into a unified theory. For temperatures around 
T ~ lOOOAvTo-the non-perturbative mu is already 
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. For a detailed discussions, we re- 



fer to H]. Large corrections to leading order mo 
have also been observed in 4D Nf — SU(2) 
simulations, using gluon mass measurements in 
the Landau gauge It remains to be seen 

whether this modification of the standard picture 
of high-temperature gauge theories has applica- 
tions in the cosmological discussion of the quark- 
hadron phase transition or in the phenomenology 
of heavy ion collisions. 
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